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Abstract 

We  consider  equilibrium  solutions  of  the  Smoluchowski  equation  for  rodlike 
nematic  polymers  with  a  Maier-Saupe  excluded  volume  potential.  The  purpose 
of  this  paper  is  to  present  a  new  and  simplified  proof  of  classical  well-known 
results:  (1)  all  equilibria  are  axisymmetric  and  (2)  modulo  rotational  symmetry, 
the  number  and  type  of  axisymmetric  equilibria  are  characterized  with  respect  to 
the  strength  of  the  excluded  volume  potential.  These  results  confirm  the  phase 
diagram  of  equilibria  obtained  previously  by  numerical  simulations  (Faraoni 
etal  1999  J.  Rheol.  43  829-43,  Forest  et al  2004  Rheol.  Acta  43  17-37,  Larson 
and  Ottinger  1991  Macromolecules  24  6270-82). 

Mathematics  Subject  Classification:  35Kxx,  70Kxx 


1.  Introduction 

The  isotropic-nematic  (I-N)  first-order  phase  transition  in  hard  rod  gases  and  liquids  is  a 
classical  topic,  which  was  first  explained  theoretically  by  Onsager  in  terms  of  an  excluded 
volume  potential  [20].  Later  Maier  and  Saupe  re-examined  the  I-N  transition  with  a 
simpler  potential  that  now  bears  their  names  [19].  The  Maier-Saupe  potential  is  a  quartic 
approximation  of  the  Onsager  potential,  which  affords  sufficient  degrees  of  freedom  to  capture 
the  hysteresis  loop  in  an  order  parameter  representation  of  isotropic  and  anisotropic  equilibria. 
Indeed,  this  idea  is  the  kernel  of  Onsager’s  insight. 
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Doi  and  Edwards  [7]  employed  a  second-moment  closure  approximation  of  the 
Smoluchowski  equation  for  the  probability  density  function  (PDF)  to  further  illustrate  the 
robustness  of  the  I-N  transition  with  coarse-grained  models  and  an  excluded-volume  potential. 
Since  then,  numerous  numerical  studies  and  semi-analytical  studies  have  provided  detailed 
bifurcation  diagrams  for  equilibrium  solutions  of  the  Smoluchowski  equation  as  well  as  its 
various  closure  approximations  [1,6, 15, 16,21]. 

These  references  yield  overwhelming  evidence  that  all  anisotropic  equilibria  of  the 
Smoluchowski  equation  are  axisymmetric.  In  [10],  the  authors  prove  this  result  for  all  second- 
moment  closure  models  with  a  quartic,  rotationally  invariant,  excluded-volume  potential.  The 
rotational  invariance  of  all  equilibria  is  the  key  ingredient  exploited  by  Onsager  to  prove 
the  I-N  transition  is  a  first-order.  This  0(3)  symmetry  of  equilibria  is  made  explicit  for 
the  Smoluchowski  equation  in  [14],  The  parametrization  of  orientational  degeneracy  of 
nematic  equilibria  is  a  prerequisite  to  rigorous  analysis  of  selection  criteria  when  full  rotational 
symmetry  is  broken  by  applied  fields  [11-13, 15].  The  proof  that  all  kinetic  model  equilibria 
are  axisymmetric  remained  open  until  quite  recently,  when  Constantin,  Kevrekidis  and  Titi, 
and  other  groups  [3-5,9, 18],  began  to  revisit  Onsager’s  seminal  papers  and  provide  modern 
rigorous  proofs  about  stationary  solutions. 

Two  groups  have  recently  proved  that  all  anisotropic  equilibrium  solutions  of  the 
Smoluchowski  equation,  or  equivalently  the  Euler-Lagrange  equation  of  the  free  energy 
for  rodlike  nematic  polymers  with  the  Maier-Saupe  excluded  volume  potential,  must  be 
axisymmetric  [8,  17],  Both  proofs  use  elaborate  estimates  on  a  scalar  integral  equation 
involving  the  eigenvalues  (or  equivalently  order  parameters)  of  the  second  moment  of  the 
PDF.  In  this  paper,  we  present  a  different  and  more  transparent  proof  for  the  same  result 
independent  of  the  work  of  [8, 17].  Our  proof  uses  elementary  calculus.  We  also  give  a  new 
proof  on  the  existence  and  the  number  of  nematic  phases  modulo  orientational  symmetry  of 
the  PDF.  Our  method  in  proving  the  number  of  nematic  phases  is  particularly  interesting  and 
may  be  extended  to  solve  many  other  problems.  It  predicts  global  properties  of  a  function  by 
studying  its  local  behaviour  at  a  hypothetical  point  with  certain  properties. 

2.  All  anisotropic  equilibria  of  the  Smoluchowski  equation  are  axisymmetric 

The  governing  equation  for  stiff  rod  nematic  polymers  interacting  by  excluded-volume  effects, 
absent  of  external  fields  and  flow,  is  the  Smoluchowski  equation  [7]: 


(1) 


where  u  is  the  unit  vector  representing  the  orientation  of  the  rodlike  nematic  polymer,  3/3 u 
the  orientational  gradient  operator  [2],  p(u,  t)  the  PDF  for  the  rodlike  nematic  polymers  with 
orientation  u  at  time  t,  D  the  rotational  diffusion  constant,  kB  the  Boltzmann  constant  and 
T  the  absolute  temperature.  In  (1),  V (u,  [p])  is  the  mean  field  excluded  volume  potential, 
which  depends  on  the  PDF  p(u,  t).  Thus,  equation  (1)  is  nonlinear  in  the  moments  of  p. 

We  adopt  the  Maier-Saupe  potential: 


F(m,  [p])  =  —bkBT(u  <g>  u)  :  (u  ®  u). 


(2) 


where  b  —  3  A/2  and  N  is  the  dimensionless  concentration,  b  can  be  viewed  as  the  strength 
of  the  excluded  volume  potential.  Here  (u  ®  u)  denotes  the  second  moment  of  the  PDF  or 
the  mean  of  the  random  variable  a®«on  the  unit  sphere  weighted  by  the  probability  density 
p(u,  t ): 


(3) 
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In  equilibrium,  the  solution  of  the  Smoluchowski  equation  is  given  by  the  Boltzmann 
distribution: 

_  exp(— V («,  [p])) 

PU  Is  expl—  y  (u,  [p]))du 

It  is  clear  that  the  equilibrium  is  completely  determined  by  the  second  moment  tensor 
(u  ®  u)  should  it  be  known.  However,  (it  <g>  it)  depends  on  the  PDF.  Moreover,  as  we  will 
show,  not  all  choices  of  (it®  it)  are  allowed.  Because  (n®  it)  is  symmetric  and  trace  one,  it  can 
be  diagonalized  by  rotating  the  Cartesian  coordinate  system  onto  its  principal  axes.  Without 
loss  of  generality,  we  assume  that  the  Cartesian  system  is  chosen  such  that  (it  ®  it)  is  diagonal. 
We  denote  it  —  (it i,  112,  u 3)  in  the  coordinate  system.  By  assumption, 


<\+n 

0 

0  ^ 

(u  <g>  u)  = 

0 

1+^2 

0 

•> 

(5) 

l  0 

0  1 

+  r3) 

where  (r\ ,  r2,  r3)  are  given 

by 

n  =  (u\)  - 

1 

3’ 

II 

"S' 

to  to 

1 

3’ 

r3  =  (u\)  -  1, 

(6) 

and  satisfy 

0  ^  |  +  r\  <  1 ,  0  ^  |  +  r2  <  1 ,  0  ^  |  +  r3  <  1 , 

r\  +  r2  +  r3  =  (u\  +  u\  +  u\)  —1  =  0. 


(7) 


The  rk,  as  eigenvalues  of  the  second  moment  of  p,  are  therefore  implicitly  defined  by  the 
PDF.  For  simplicity,  we  absorb  k^T  into  b.  The  Maier-Saupe  potential  can  be  recast  in  these 
coordinates, 

y  (it)  —  —b(r\u\  +  r2il  +  >"31/3  +  |). 

The  equilibrium  PDF  is  expressed  as  follows: 

,(«)  =  exp  [fclfl„|+.2.J|  +  ,-,„|>]  (8) 

Js  exp[h(rinf  +  r2«2  +  #"31/3)]  dt/ 

In  the  analysis  below,  we  treat  b  as  a  variable  (it  is  either  the  normalized  volume  fraction 
of  rods  or  the  reciprocal  of  temperature)  and  view  the  PDF  and  all  averages  (expectations)  as 
functions  of  b.  Specifically,  we  define  functions  as  follows: 

Rm  =  (u\)-  i,  R2(b)  =  (u\)  -  R3(b)  =  (uj)  - 

where  the  average  is  taken  with  respect  to  the  PDF  given  in  (8).  If  (n,  >'2,  '3 )  corresponds  to 
an  equilibrium  for  b  =  ho,  then  functions  A’|  (/;).  Ri(b)  and  R3(b)  must  satisfy 


/MO)  =  0,  tf2(0)  =  0,  R3(  0)  =  0, 

Ri(b0)  =  n,  R2(b0)  =  r2,  R3(b0)  =  r3. 


It  is  important  to  note  that  Rk(b)  is  a  function  of  b  while  rk  is  a  fixed  value. 

To  simplify  the  presentation,  we  introduce  some  shorthand  notation: 

hi  —  Mj  —  (Mj),  h2  =  ul  ~~  (“2)’  ^3  =  u\  —  (u\).  (10) 

We  note  that  h  1 ,  A  2  ■  h3  are  random  variables  and  functions  of  b.  They  satisfy  the  following 
constraints: 


hi  +  h  2  +  /z  3  =  0 


and 


(hi)  =  (h2)  =  (h3)  =  0. 


(ID 
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Any  solution  of  (1)  with  three  distinct  eigenvalues  of  the  second  moment  is  called  biaxial. 
PDFs  with  two  distinct  eigenvalues  of  the  second  moment  are  called  uniaxial  (axisymmetric). 
The  isotropic  state  corresponds  to  r\  =  r2  =  r3  =  0. 

Theorem  1.  All  anisotropic  equilibria  of  the  Smoluchowski  equation  with  the  Maier—Saupe 
excluded  volume  potential  are  axisymmetric.  The  isotropic  state  is  a  non-degenerate 
equilibrium  for  all  non-zero  strengths  of  the  potential. 

To  prove  this  theorem,  we  need  a  lemma. 

Lemma  1.  Suppose  r\  >  r2.  Then  for  b  >  0,  we  have 

(h3(h,  -hf))  <0, 

where  h\.  I12,  h3  are  defined  in  (10). 

Remark  1.  (h3h\)  =  (( u\  —  { u\))(u\  —  { u\ >)>  is  the  correlation  of  w2  and  u\. 

Remark  2.  This  lemma  does  not  impose  any  condition  on  r3. 

Remark  3.  We  can  permute  (n,  r2,  r3)  and  apply  this  lemma  to  any  pair  of  (r\,  r2,  r3). 

For  example,  if  r2  >  r3,  then  we  have 
(hl(h2-h3))<  0. 


Proof  of  lemma  1.  To  prove  this  lemma,  we  first  point  out  the  fact  that  the  function 
/02?r  cos  20  exp(r  cos  20)  d 0 


g(r )  = 


fQ  n  exp(r  cos  20)  d 0 


(12) 


is  an  increasing  function  of  r,  which  was  proved  in  [4].  We  rewrite  (h3(h\  —  h2))  as 
(h3(hi  -  h2))  =  (( u\  -  (u\))(u\  -u\  +  (iq)  -  (u\))) 

=  ((«3  -  («3»(«1  -  ul))- 

We  select  the  axis  associated  with  u3  as  the  z-axis  and  establish  a  spherical  coordinate 
system  (cp,  0)  where  f  is  the  polar  angle  and  0  is  the  azimuthal  angle.  We  express  everything 
in  terms  of  (</>,  0): 

u\  —  u\  =  sin2  4>  cos2  0  —  sin2  4>  sin2  0  =  sin2  cf>  cos  20, 
ul  —  cos2  f. 

The  Maier-Saupe  potential  is 

V(<p,0)  —  —b(r\  sin2  f  cos2  0  +  r2  sin2  (p  sin2  0  +r2  cos2  <p) 


=  —b(a  1  sin2  <p  +  a2  sin2  cp  cos  20  +  r3  cos2  cp), 


where 


a  1  = 


r  1  +  r2 


The  probability  density  is 
p(cp,0)  = 


r\  ~r2 

a2  =  — - —  >  0. 


exp[f>(«i  sin2  cp  +  a2  sin2  (p  cos  20  +  r3  cos2  cp)] 


fa  fo  ^  exP[^(fli  sin2  cp  +  a2  sin2  (p  cos  20  +  r3  cos2  <p)\  d 0  sin  (p  d (p 
Notice  that  the  probability  density  satisfies 


pin 

Jo  ' 


cos  2 0p(<p,  0)  dO  —  g(ba2  sin“  cp) 


p2n 

JO 


p(cp,  0)  d 0, 
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where  function  g(r)  defined  in  (12)  is  an  increasing  function  of  r.  Because  b  >  0  and  a2  >  0, 
g(bci2  sin2  0 )  is  an  increasing  function  of  sin2  0. 

Let  cos2  0 o  =  (u\)  =  (cos2  0).  Writing  the  average  in  terms  of  ( 0 ,  0),  we  obtain 

((«3  -  (m3»(“i  -  ul)) 


f 


=  I  (cos  0  —  cos  <po)  sin-  0 

/ 0 


L  ' 


r 

Jo 


=  I  (cos2  0  —  cos2  </>o)  sin2  </)g(fifl2  sin2  0) 


cos  2 0p((p,  9)  d 6  )  sin  0  d0 

2  TT 

0 


p(0,  9)  d&  )  sin0  d0 


=  ((cos2  0  —  cos2  0o)  sin2  <pg(ba2  sin2  0)) 

"7  9  9  9  9 

=  ((cos  0  —  cos  0o){sin  <pg(ba2  sin-  0)  —  sin"  0og(fia2  sin  0o)})  <  0. 

Here  we  used  that  g(ba2  sin2  0)  is  an  increasing  function  of  sin2  0  and  cos2  0=1  —  sin2  0 


is  a  decreasing  function  of  sin-  0. 


□ 


Proof  of  theorem  1.  The  PDF  (8)  depends  on  b.  Its  derivative  with  respect  to  b  is  given  by 

-ttP(m)  =  [(nuj  +  r2t<2  +  riu\)  -  {riu\  +  r2u\  +  r3M2)]  •  p(u) 
a  b 

=  {r\h i  +r2h2  +  r3h3)  ■  p(u). 

The  derivative  of  Rk  (b)  with  respect  to  b  is  given  by 

R'k(b)  =  (u\(r  \h\  +  r2h2  +  r3h3))  =  {hk(nhi  +  r2h2  +  r3h3)) , 

taking  into  account  the  fact  (hi)  =  ( h2 )  =  (h3)  =  0. 

We  prove  the  theorem  by  contradiction.  Suppose  there  is  a  set  of  distinct  (n,  r2,  r3) 
satisfying  equation  (6)  with  probability  density  (8)  for  b  =  bo  >  0.  Without  loss  of  generality, 
we  assume  r\  >  r2  >  r3.  Because  r i  +  r2  +  r3  =  0,  we  have  r\  >  0  and  r3  <  0. 

We  consider  the  function 

F(b)  =  r2Ri(b)  -  r1R2(b). 

Equation  (9)  implies  that  F( 0)  =  F(bo)  =  0.  To  establish  the  contradiction,  we  are  going 
to  prove  that 

F'(b)>  0  for  b  >  0. 

Taking  the  derivative  of  F(b)  with  respect  to  b,  we  have 
F'(b)  =  r2R[(b)-riR'2(b) 

=  r2(hi(rihi  +  r2h2  +  r3b3))  -  ri(h2(r\h\  +  r2h2  +r3h3)). 

Using  hi  =  —(h2  +  h3)  and  h2  =  —  (hi  +  h3),  we  obtain 
F'(b)  =  r2(hi[(r2  -  n)h2  +  (r3  -  n)h3])  -  ri(h2[(ri  -  r2)hi  +  (r3  -  r2)h3]) 

=  r2(r2  -  ri)(hih2)  +  r2(r3  -  ri)(hxh3)  -  ri(ri  -  r2)(hih2)  -  n(r3  -  r2)(h2h3). 

Splitting  the  term  of  (hih3)  into  two  using  the  identity  (r3  —  ri )  =  (r3  —  r2)  +  (r2  —  r i), 
combining  two  terms  of  (hih2)  and  using  r3+r2  =  —r3,  we  arrive  at 

F'(b)  =  r2(r3  -  r2)(hih3)  +  r2(r2  -  ri)(h\h3)  -  r3(r2  -  ri)(hih2)  -  ri(r3  -  r2)(h2h3) 

=  (r3  -  r2)[r2(hih3)  -  ri(h2h3)]  +  (r2  -  ri)[r2(hih3)  -  r3(hih2)]. 

Applying  the  identity 

ctiPi  -  c(2i 62  =  j[(q:i  -  a2)(Pi  +  £2)  +  («i  +  u2)(Pi  -  P2)] 
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and  using  r\  +  r2  =  —  r3  and  r2  +  r3  =  —  r\,  we  finally  obtain 
F'(b)  =  i{(r3  -  r2)[(r2  -  n)(h3(h i  +  h2))  +  (r2  +  n)(h3(hx  -  /t2)>] 

+(7*2  -  n)[(r2  -  r3){h\{h3  +h2 ))  +  (r2+r3)(h1(h3  -  h2))]\ 

=  k{(r3  -  r2)(r2  -  n)\{h3(hi  +  h2))  -  (hl(h3+h2))] 

Hr 3  -  f*2 )(  r3)(h3(hi  -  h2 ))  +  (r2  -  7-i)(-n)(/!i(/i3  -  ^2))} 

=  3 {(^1  -  r2)(r2  -  r3){h2(h3  -  hi))  +  (r2  -  r3)(-r3){h3(h2  -  /zi)> 

+  (7*1  -  7*2)7*l(/7l(/73  -  /r3))}- 

In  the  above,  all  coefficients  are  positive,  and  lemma  1  dictates  that  all  correlation  terms 
are  positive  as  well.  Therefore,  we  conclude  that 

F'(b)>  0  for  b  >  0, 

which  clearly  contradicts  (9).  Therefore,  there  cannot  exist  three  distinct  r k  satisfying  equation 
(6)  with  PDF  (8).  That  is,  the  equilibrium  solutions  of  the  Smoluchowski  equation  are 
parametrized  by  three  r k  of  which  at  least  two  must  be  equal. 


3.  Number  of  equilibria 


Given  that  all  equilibria  are  axisymmetric,  we  now  investigate  their  multiplicity  (modulo 
rotational  invariance).  We  select  the  axis  of  symmetry  as  the  z-axis.  We  next  establish  a 
corresponding  spherical  coordinate  system.  Let  (p  be  the  polar  angle  measured  from  the  z-axis, 
0  ^  (p  ^  jt.  Let  0  be  the  azimuthal  angle,  the  angle  on  the. ry-plane  measured  counterclockwise 
from  the  x-axis,  0  ^  6  ^  2jt.  We  express  u  in  terms  of  ( (p ,  6): 

u=  (mi,  u2,  u3)  —  (sin(0)  cos(0),  sin(</))  sin(0),  cos(</>)).  (13) 


Because  of  axisymmetry,  (5)  becomes 

/  5  (sin2  (</>)}  0 

0  i  (sin2  (0)) 


V 


0  \ 

0  ,  (14) 

0  0  (cos  2((p))J 

and  the  Maier-Saupe  potential  simplifies  accordingly: 

(u  0  u)  :  ( u  ®m)  =  \  sin2 (0)  (sin2 (0))  +  cos2 {(p)  (cos2 ((p)) 

=  |(3(cos 2(0)>  -  1)  cos2 (</>)  -  ^ ((cos2 {(p))  -  5)  +  5.  (15) 

Notice  that  the  last  two  terms  on  the  right  side  of  (15)  are  constants,  independent  of  (p. 
Thus,  the  Maier-Saupe  potential  (up  to  a  constant)  is 

V(<p)  =  -Z?i(3(cos2(c/>)}  -  1) cos 2W>).  (16) 

Here  we  introduce  a  new  variable  following  Constantin  and  Vukadinovic  [5] : 

r  =  Z?i(3(cos2(0)>  -  1),  (17) 

where  we  note  that  i(3(cos  2((p))  —  1)  is  the  Flory  order  parameter. 

We  rewrite  the  potential  as 

V {(p)  =  —r  cos2 (0),  (18) 


so  that  the  Boltzmann  distribution  (4)  in  spherical  coordinates  has  the  form: 

exp(r  cos2  (</>)) 

2i x  exp(r  cos2 (</>))  sin{(p)  dcp 


P((p)  = 


(19) 
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We  emphasize  that  r  has  to  satisfy  equation  (17)  with  PDF  given  by  (19): 

1  / 3  f*  cos 2(0)  exp(r  cos2  ((/>))  sin(</>)  d (f> 

2y  exp(r  cos2 ((/>))  sin (</>)  d0 

Using  w  —  cos(0),  (20)  becomes 

1  f}  (3w2  —  1)  exp(rw2)  dw 

r  —  b - j - . 

2  f.  exp (rw2)  dw 

Integrating  by  parts  and  collecting  terms  to  the  left,  we  obtain 

r(l  -  bf{r))  =  0,  (21) 

where  f(r)  is  given  by 

def  L  w2(l  —  W2)  exp(rw2)  dw 

fir)  =  ^ ^ - ’—f. - .  (22) 

J0  cxp(rw-)  dw 

The  equilibrium  equation  (21)  is  a  nonlinear  integral  equation  for  the  scaled  Flory  order 
parameter  r.  For  each  solution  of  equation  (21),  the  corresponding  equilibrium  solution  of 
equation  ( 1 )  is  given  by  ( 1 9).  Therefore,  the  number  of  equilibrium  solutions  equals  the  number 
of  solutions  of  (21).  We  first  note  that  r  —  0  solves  (21)  for  all  values  of  b.  When  r  =  0,  from 
(18)  it  follows  that  r\  —  ri  =  =  0  and 

(m  ®  u)  =  j/, 

which  translates  back  to  the  PDF  as  the  isotropic  solution  p  —  1  /Ait . 

For  b  =  0,  r  =  0  is  the  only  solution  of  (21).  For  b  >  0,  non-zero  solutions  of  (21)  must 
satisfy 

l  =  fir).  (23) 

b 

Theorem  2.  f  (r)  defined  in  (22)  has  the  following  properties. 

(1)  /( 0)  = 

(2)  0  <  fir)  <  f 

(3)  limr_5._00  f(r)  =  0  and  limr^+00  fir)  =  0. 

(4)  There  exists  r *  >  0,  such  that  fir*)  —  0,  fir)  >  0  for  r  <  r*  and  fir)  <  Ofor  r  >  r*. 

Proof  of  theorem  2.  To  prove  this  theorem,  we  first  point  out  a  fact  of  calculus  of  one  variable. 

Suppose  fir)  satisfies  the  property  that  f'(rfi)  —  0  implies  /"(ro)  <  0.  If  fir)  attains  a 
maximum  at  r* ,  then  we  have 

fir)  >0  for  r  <  r*  and  fir)  <  0  for  r  >  r*. 

As  we  will  see,  this  simple  fact  of  calculus  plays  a  very  important  role  in  the  proof  of 
property  (4)  below.  It  enables  us  to  predict  the  global  behaviour  of  a  function  by  studying 
its  local  behaviour  at  a  hypothetical  point  with  certain  properties.  Now  we  prove  the  four 
properties  listed  one  by  one. 

(1)  /( 0)  =  /(,'  w2(l  -  w2)  dw  =  2  -  i  =  2L. 

(2)  0  <  w2(l  -  w2)  <  \  for  u>  €  (0,  1)\  { -4 }  lead  to 

(l  w2(  1  —  w2)  exp(rw2)  dw  1 

0  <  fir)  =  2° - - L -  < 

fQ  exp(rw2)  dw  4 
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(3)  To  prove  limr_  ,  -oc  fir)  —  0,  we  consider  for  any  r  <  0  and  1  >  s  >  0, 

/  exp(rw2)  dw  —  I  exp[r(e  —  s)2]  d.s  =  /  expire2  —  2res  +  rs 2]  ds 

Jo  Jo  Jo 

^  exp  (re2)  -  /  exp(— rs2)  d.v  ^  /  exp(rw2)dw-  /  exp(— rs2)  d.v. 

Jo  Je  Jo 

Here  we  have  used  —2 res  +  rs2  ^  —2 rs2  +  rs2  —  —rs2. 

Since  limr  _>_<*>  So  exP(  —rs2)  d.v  =  oo,  we  obtain  the  result  below. 

For  any  e  >  0,  there  exists  an  M  such  that  r  <  —M  implies 


f 


exp(rur)  d  w  ^  - 
£ 


L.r 

£  Je 


exp(rw“)  dw. 


It  follows  that 


0  <  fir)  < 


e2  /0£  exp(rw2)  dw  +  (1/4)  /f'  exp (rw2)  dw 
/0e  exp(rw2)  dw 


<  e 


Thus,  we  have  limr^_oo  /(r)  =  0. 
lim,.^+0O  f{r)  —  0  is  proved  in  a  similar  way. 

(4)  Because  linv^ioo  fir)  —  0  and  fir)  >  0,  fir)  attains  a  maximum  at  r*.  Using  the 
simple  fact  of  calculus  we  pointed  out  above,  we  see  that,  to  prove  property  (4),  we  only  need 
to  prove  that  f'irf)  —  0  implies  /"(ro)  <  0. 

For  the  simplicity  of  presentation,  we  consider  random  variable  IT  with  probability  density 

exp(rw2) 

pico)  =  - . 

j0  exp(rw2)dw 

Then  we  have  fir)  =  ( IT2(1  —  W2)).  The  derivative  of  p(w)  with  respect  to  r  is 


J-piw)  =  iw2  -  { W 2»  •  piw). 
dr 

Taking  derivatives  of  /(r)  =  {W2{  1  —  W2)),  we  have 

fir)  =  {W2i\  -  IT2)(IT2  -  (W2))) 

=  (W2i  1  -  W2)W2)  -  (W2i  1  -  W2))  ■  {W2) 

=  ({IT2(1  -  W2)  -  (w2i  1  -  VU2)>}VU2> 

=  ({(W2(l  -  IT2))  -  w2il  -  W2)}(1  -  IT2)), 
fir)  =  ({(IT2(1  -  W2))  -  W2i  1  -  W2)}(1  -  W2)iW2  -  (IT2))) 

+  (T(it2(i-it2)))-((1-it2)> 

=  «(IT2(1  -  IT2))  -  w2il  -  W2)}IT2(1  -  W2)) 

-  ({(IT2(1  -  W2))  -  w2i  1  -  w2)}i  1  -  W2))  ■  (IT2)  +  fir)  •  (1  -  (IT2)) 
=  (W2(l  -  IT2))2  -  ({ir2(l  -  IT2)}2)  +  fir)  ■  (1  -  2(IT2)). 

Here  we  have  used 

J-(W2H-W2))  =  fir) 
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r 

Figure  1.  Graph  of  function  f(r)  defined  in  (22). 


and 

({(W2(l  -  W 2)>  -  W2(  1  -  W2)}(  1  -  W2))  = 

Suppose  f'(ro)  =  0.  Then  we  have 

f"(r0)  =  (W2(  1  -  W2))2  -  ({W2(\  -  W2)}2)  =  -var{W2(l  -  W2)}  <  0. 


Then  it  follows  that  f'(r)  >  0  for  r  <  r*  and  fir)  <  0  for  r  >  r* . 
Finally  we  have 


which  implies  r*  >  0. 


r*  and  f(r*)  can  be  estimated  numerically  using  a  nonlinear  equation  solver  (such  as  the 
bisection  method)  and  a  numerical  integration  method  (such  as  Simpson’s  method).  Our 
numerical  results  show 

r*  «  2.178  287  9748,  f(r*)  »  0.148  555  599  922  54. 

A  plot  of  f(r)  is  shown  in  figure  1.  □ 

Corollary.  Let  b*  =  1  /f(r*)  &  6.731486  3965. 

(1)  For  b  <  b* ,  (21)  has  one  solution,  r  —  0. 

(2)  For  b  —  b*,  (21)  has  two  solutions,  r  =  0  and  r  =  r*. 

(3)  For  b*  <  b  <  4^,  (21 )  has  three  solutions  r  —  0,  r  =  r\  and  r  =  r2  where  0  <  r\  <  r* 
and  ri  >  r*. 

(4)  For  b  —  (21)  has  two  solutions  r  =  0  and  r  —  ri  where  r2  >  r* . 

(5)  For  b  >  y,  (21)  has  three  solutions  r  =  0,  r  =  r\  and  r  =  r2  where  r\  <  0  and  >  r*. 

Remark.  Our  numerical  value  of  b*  is  slightly  different  from  the  numerical  value  of  6.73 1  393 
reported  in  [17]. 
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4.  The  2D  case 


Our  approach  also  works  for  the  2D  Smoluchowski  equation.  In  the  proof  given  by  Constantin 
and  Vukadinovic  [5]  for  the  2D  problem,  the  key  step  is  to  prove  that  the  function 


/0?7r  cos"  4>  exp(r  cos  <p)  &<p 


f(r)  = 


/02?r  exp(r  cos  <p)  &<p 


is  an  increasing  function  of  r.  Constantin  and  Vukadinovic  [5]  proved  that  fir)  is  strictly 
increasing  by  viewing  it  as  the  solution  of  a  differential  equation.  Here  we  use  the  simple  fact 
of  calculus  pointed  out  in  the  proof  of  theorem  2  above.  Again,  the  main  idea  of  our  approach 
is  to  predict  the  global  behaviour  of  a  function  by  studying  its  local  behaviour. 

We  first  prove  that  fir)  satisfies  the  property  that  /'(ro)  =  0  implies  /"(ro)  >  0. 

We  consider  random  variable  <t>  with  probability  density 

...  exp(r  cos/) 

Pi<t>)  =  - ■ 

L  exp(r  cos/)  d/ 

Then  we  have  fir)  =  (cos2  <t>).  The  derivative  of  p((p)  with  respect  to  r  is 
d 

— p(/)  =  (cos  /  -  (cos  0>»  ■  p  if). 
dr 


Taking  derivatives  of  fir)  =  (cos2  <J>),  we  have 


fir) 

f"(r) 


(cos2  0(cos  <l>  —  (cos  O})} 

((cos2  $  —  (cos2  <J>})  cos  <P), 

((cos2  —  (cos2  $})  cos  cJ>(cos  $  —  (cos  $})}  —  ^^-(cos2 
((cos2  $  —  (cos2  <J>})  cos2  O) 

—  ((cos2  <J>  —  (cos2  <J>>)  cos  <p)  •  (cos  <h>  —  f'(r)  ■  (cos  0> 
(cos4  <f>)  —  (cos2  O)2  —  2/'(r)  •  (cos  O). 


(cos  cp) 


Here  we  have  used 


d 

dr 


(cos2  O)  =  f(r) 


and 


((cos2  <1>  —  (cos2  <h})  cos  $)  =  /'(r). 

Suppose  /'(ro)  =  0.  Then  we  have 

/"(r0)  =  (cos4  <5>)  —  (cos2  T>)2  =  var  {cos2  O}  >  0. 

Now  we  show  that  /'(r)  >  0  for  r  >  0. 

Because  /'( 0)  =  0  we  obtain  /"( 0)  >  0  ,  which  means  that  /'(r)  >  0  for  r  e  (0,  e). 

Let  ro  =  sup{^|g  >  0  and  fir)  >  0  in  (0,  q)}.  If  ro  =  oo,  then  it  is  true  that  fir)  >  0 
forr  >  0.  If  ro  <  oo,  then  we  have  /'(ro)  =  0  and  consequently  /"(ro)  >  0.  That  means 
fir)  <  0  for  r  e  (ro  —  s,  ro),  which  contradicts  the  definition  of  ro- 
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5.  Conclusion 

We  have  presented  three  new  proofs  for  the  type  and  number  of  the  equilibrium  solutions 
of  the  Smoluchowski  equation  with  Maier-Saupe  excluded  volume  potential  in  2D  and  3D, 
respectively.  Although  the  conclusions  are  not  new,  the  approach  is  new  and  systematic  and 
has  potential  applicability  to  other  questions  such  as  selection  criteria  for  stationary  solutions 
of  a  Smoluchowski  equation  in  the  presence  of  axisymmetric  external  flow  and/or  external 
fields  coupled  to  an  excluded-volume  potential. 
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